Comment on "New Methods for Old Coulomb Few-Body Problems 



V. S. Zotev 

Los Alamos National Laboratory, Biophysics Group, MS D454, 



Los Alamos, NM 87545, USA 



in 
o 
o 

c 
3 



43 

Oh! 



o : 
^— > 

O ■ 

43 ' 

m : 
> ■ 



m ■ 

o ■ 
in ; 

o . 

o : 
• i— i 

43 : 

Ph. 



X 



In this Comment on the above mentioned paper by F. E. Harris, A. M. Frolov, and V. H. Smith, 
we briefly review our contributions to development of new methods for solution of the Coulomb four- 
body problem. We show that our research group, headed by Prof. T. K. Rebane, had a priority in 
using the fully correlated exponential basis for variational calculations of four-body systems. We also 
draw attention to the fact that our group subsequently implemented a more advanced method, which 
uses highly efficient exponential-trigonometric basis functions for solution of the same problem. 



The paper "New Methods for Old Coulomb Few-Body 
Problems" pj deals with analytical and methodological 
issues involved in the implementation of two methods - 
the Hyllcraas method and what the authors call the "ex- 
ponential ansatz" approach. It is stated in the abstract 
that the second method "has only recently been imple- 
mented for four-body systems" [jj. The last section of 
the paper is entitled "Preliminary Results" p|. 

The main purpose of this Comment is to call attention 
to our works in this field and show that the method using 
a large fully correlated exponential basis for solution of 
the Coulomb four-body problem was in fact implemented 
and successfully applied by our group ten years ago. 

A general method for analytic evaluation of four- 
particle integrals was developed by Fromm and Hill in 
their well-known pioneering work |2( in 1987. This 
method makes it possible to perform direct variational 
calculations of four-body systems using exponential ba- 
sis functions of all six interparticle separations. The first 
such calculation of two model systems was carried out 
by Rebane and Yusupov in 1992 Q. Rebane showed in 
1993 how the number of different types of four-body in- 
tegrals that need to be computed can be reduced from 43 
to only seven Q . The first variational calculation of the 
positronium molecule e + e~e + e - in the fully correlated 
exponential basis was performed by Rebane, Zotev, and 
Yusupov in 1996 [{|. It employed 90 symmetrized expo- 
nential functions. This work was followed by calculations 
of eight symmetric mesomolecules in the basis of 240 ex- 
ponentials 0- In 1998, the original method by Fromm 
and Hill was generalized by Zotev and Rebane to allow 
evaluation of four-particle integrals with complex expo- 
nential parameters. This enabled us to perform the first 
variational calculation of the positronium molecule and 
three other systems using the basis of extremely flexi- 
ble exponential-trigonometric functions p| • It was shown 
that a single exponential-trigonometric function is as effi- 
cient as seven exponentials in the case of e + e~e + e~ , and 
is flexible enough to allow direct variational computation 
of even adiabatic four-particle systems such as the hy- 
drogen molecule Q . The exponential-trigonometric basis 
was initially introduced by our group in the three-body 
problem y|. A detailed description of our methods us- 
ing the exponential and exponential-trigonometric basis 
functions for direct variational solution of the Coulomb 
four-body problem was published in 2002 . 



Harris, Frolov, and Smith implemented the exponen- 
tial basis and reported a calculation of the positronium 
molecule with two exponential functions in 2003 |10| . 
They used a modified Fromm and Hill formula |ll|. and 
argued that the advantage of their approach was "better 
formal treatment of termwise singularities" [T(il |. They 
wrote in the discussed paper that their refinements en- 
abled them to find "a two-configuration wavefunction 
with energy -0.513237 hartree, comparable to the result 
obtained by Rebane et al. using four configurations" p|. 

Our extensive practical experience shows that the 
problem of termwise singularities is not serious, because 
parameters of basis functions can easily be chosen to 
avoid singularities in all the integrals to be computed. 
Our calculations of the positronium molecule [5j and me- 
somolecules |6j with the proper parameter selection ex- 
hibited perfect numeric stability and steady convergence. 

The authors' comparison of their two-function calcula- 
tion and our four- function result is performed incorrectly. 
In their calculation [Toj , the two functions were optimized 
simultaneously. In our work p|, the basis was expanded 
gradually, and basis functions were optimized one at a 
time. Naturally, the simultaneous optimization provides 
better results. Moreover, our 1996 work included not just 
four, but a total of 90 exponential functions, and yielded 
a variational energy value of -0.515920 hartree 

Finally, we have to mention the following. Frolov and 
Smith pj, Frolov p]|, Frolov and Bailey pj, when dis- 
cussing the Fromm and Hill four-body integrals in their 
2001-03 articles, did not reference our works at all. Sim- 
ilarly, Frolov made no reference to our original works 
[8| in his 1998-03 papers, which used the exponential- 
trigonometric basis in the three-body problem [l5j . In 
the discussed paper Harris, Frolov and Smith cite 
our 1996 calculation of the positronium moleculeJaK but 
do not reference any subsequent publications |fj, 0, 
particularly the 2002 article [j| , which provides a detailed 
review of our methods. Frolov and Smith follow the same 
line in another recent paper p"fj| . The most advanced and 
accurate method, which uses highly flexible exponential- 
trigonometric functions for calculations of four-particle 
systems with arbitrary particle masses 0, E3 , has never 
been mentioned by Harris, Frolov, and Smith. Even the 
latest and most comprehensive review article by Harris 
[l7T | makes no explicit mention of this method. We leave 
these facts without comment. 
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